Universality and tree structure of high energy QCD 
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Using non-trivial mathematical properties of a class of nonlinear evolution equations, we obtain 
the universal terms in the asymptotic expansion in rapidity of the saturation scale and of the 
unintegrated gluon density from the Balitskii-Kovchegov equation. These terms are independent 
of the initial conditions and of the details of the equation. The last subasymptotic terms are new 
results and complete the list of all possible universal contributions. Universality is interpreted in a 
general qualitative picture of high energy scattering, in which a scattering process corresponds to a 
tree structure probed by a given source. 



1. Introduction 

The energy dependence of hard hadronic cross sections 
is encoded in the evolution of the parton distributions 
with the rapidity Y. If the relevant scale Q is large 
enough, like in the case of deep-inelastic scattering, for 
which Q is the virtuality of the photon, the latter can 
be computed in perturbative QCD. The Balitskii-Fadin- 
Kuraev-Lipatov (BFKL) equation resums the leading 
Fcynman diagrams, in the form of a linear integrodiffer- 
ential evolution equation for the gluon density. Its solu- 
tion exhibits an exponential growth with Y. At larger 
Y, maybe already attainable at present colliders the 
BFKL equation breaks down because finite parton den- 
sity or saturation effects become sizeable 3, 4]. It is 
replaced by nonlinear evolution equations which aim at 
describing such effects 0, 0, • The appropriate equa- 
tions depend on the physical observables and on the level 
of approximation of QCD. They could also include non- 
perturbative contributions 8]. Our aim is to show that 
minimal assumptions on the structure of saturation equa- 
tions lead to quite general, i.e. "universal" , properties of 
the solutions. As a case study, we will consider the Bal- 
itskii-Kovchegov (BK) equation 0, 

In the large Nc limit, at fixed QCD coupling a = 
asNc/n, and for a homogeneous nuclear target, the mea- 
sured gluon density Af{k, Y) at transverse momentum k 
obeys the BK equation, which reads 



(1) 



Ol stands for the partial derivative with respect to log fc^ 
and the function x is the BFKL kernel xil) ~ 21/1(1) — 
■0(7) — '0(1 — 7). The nonlinear term in the r.h.s. corre- 
sponds to the saturation term. For its relative simplicity. 
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the BK equation is the prototype of saturation equa- 
tions. 

However, besides the mathematical challenge posed 
by the nonlinear character of saturation equations, 
there is up to now no fully consistent way to include 
saturation corrections in the perturbative expansion of 
QCD. Furthermore, initial conditions are required at 
some low rapidity Yq, for all values of A:: it is a priori 
not clear how dependent on them the solution of Eq. 
is, in particular for small k ~ Aqcq, which lies beyond 
the known perturbative domain. Such sensitivity on 
the infrared has already been a major drawback to the 
BFKL equation 0. One purpose of this Letter is to 
explain in what respect the solutions to the BK equation 
are independent of the details of both the nonlinear 
evolution equation and of the initial conditions. Let us 
call this property universality. 

2. Universal properties of the solutions of the BK 
equation 



We are going to use general methods for solving non- 
linear evolution equations [lol lll| . The mathematical 
results that will be relevant for the BK equation require 
that three main conditions be fulfilled [ij. (i) = 
should be a linearly unstable fixed point, i.e. any small 
fluctuation of M grows to infinity from the linearized evo- 
lution equation; (ii) the effect of the nonlinearity is to 
damp this growth when Af ^ 0{1); (Hi) the initial con- 
dition must be "sufficiently" steep at large k. Properties 
(i) and (ii) stem from the parton model in general and 
are verified by the BK equation in particular. In the 
case of the BK equation, as noticed in Refs. [l^ Il3| . 
the constraint (Hi) means A/'(fc,lo) *C l/k'^''", with jc = 
0.6275 • • • . Color transparency of perturbative QCD im- 
plies Af{k,Yo) ^ 1/k'^ at large k, which fulfills this con- 
dition. 

The class of equations defined by the latter constraints 
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admit traveling wave solutions""^ , i.e. at large Y, Af is a 
uniformly translating front: 



AA(fc, y ) ^ AAoo ( log fc2 _ log Ql{Y)) 



(2) 



The mathematical treatment of Eq. consists in con- 
sidering asymptotic expansions in Y of both the satura- 
tion scale (related to the translation of the traveling wave 
front with rapidity Etnd of the gluon distribution A/" 
(which is the front profile It leads to a hierarchical 

system of nested ordinary differential equations. At each 
order of the expansion, the integration constants are de- 
termined by a matching procedure \vi\ . The correspond- 
ing general equations can be straightforwardly translated 
to the case of the BK equation. 

Let us first discuss the saturation scale Qs{Y). In 
terms of the parameters of the BK equation, the sys- 
tematic expansion^ of Qs{Y) in Y reads 



logQliY)^a^Y 



Ic 



T^iogy 

27c 



2ti 



1 



Ic V ax" (7c) \/r 



+ 0{1/Y). (3) 



7c is the solution of the implicit equation x(.lc) = 
7cx'(7c), see Ref. Q. The successive orders contribut- 
ing to dlogQ^giY) / dY are displayed in Fig. 1. 

A few comments are in order. The first term in Eq. ^ 
is the dominant evolution term [l6|. The second term 
[itI I reflects the delay in the formation of the front [l3j |. 
The third term is a new result for the solution of the 
BK equation. These three terms are universal. They 
are the only possible universal terms in this asymptotic 
expansion. Indeed, one sees that by performing a shift 
in the rapidity Y Y -\-Yq m. Eq. that amounts to 
modifying the initial conditions, one gets additional 1/Y 
terms, while the lowest order terms are not affected by 
the shift. 

We now turn to the discussion of the properties of the 
front. The whole method to get solutions relies on the 
matching between two different expansions of the wave 
front above the saturation scale: the so-called "front 
interior" and "leading edge" expansions the 
"front interior" , the relevant small parameter is z = 
\og[k'^ /Ql{Y))/ ^2ax"{lc)Y < 1, i.e. one stands at 
transverse momenta k near the saturation scale QgiY). 



^ Rigorous mathematical results are ava ilab le for the Fisher 
and Kolmogorov-Pctrovsky-Piscounov Il5l equat ion (related to 
the BK equation in the diffusive app rox imation |12| '). and have 
been extended to more general cases 

^ Note that the scale of Qs{Y) is included in the definition. Indeed, 
the saturation scale Qs{Y) can be defined, e.g., in such a way 
that J\f{Qa(Y), Y) = A/q, where A/q is a given constant. QsiY) 
is the position of the traveling wave front. All universal terms 
are independent of A/q . 
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FIG. 1: The effective dependence d\ogQl{Y)/dY of the sat- 
uration scale upon the rapidity Y. a was set to 0.2. The 
truncations of the asymptotic expansion Eq. Q to order 1, 
l/Y and l/Y''^^^ are shown. 



J\f is expanded about its asymptotic shape A/'oo, see 
Eq. 10). In that region, Af ^ I and the full nonlinear 
equation must be considered. In the "leading edge" de- 
fined by z ~ 1, the wave front is expanded in powers of 
1/Vy : it is the transition region, where the front forms. 
Going from the front interior to the leadi ng e dge, the 
asymptotic front profile Moo crosses over to 



Af{k,Y) = C I ^' .-z^ ) i.„r,.2 



e-^ X |7,log[fcVQ^(y)] 

X"{lc) 



h C2 + (^3 - 2C2 

27cX^%c) 1 r 2] 

6V^{1-,F^[~i,I;z'])z + 0{1/Vy)\ (4) 



where Ci and C2 are two integration constants and iFi , 
2F2 are hypergeometric functions. 

The first term can be identified to A/'oo x near 
the front, and reproduces the known results [l^ ll6lll^ : 

A/'oo induces so-called geometric scaling [l8|, while 

2 

e^^ describes the diffusive formation of the front, and 
induces a specific pattern of geometric scaling violations. 
All other terms appearing in the expansion Q are 

new terms that represent the 0{1/VY) corrections to 

2 

A/00 X . Note that they depend also on the third 
derivative X^^Hlc), at variance with the new term of 
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order 1 / \/Y in the expression for the saturation scale Q . 




FIG. 2: The reduced front profile M{k,Y) x {k^ /Ql{Y)y-= 
in the leading edge . The curves are for y = 10 (lower 
curves) to F = 50 (upper curves), ct — 0.2 and the con- 
stants in Eq. ||1J were set to Ci = 1, C2 ~ 0. The lead- 
ing edge expansion I^J is truncated at level 0{VY) (dashed 
lines) and 0(1) (solid lines). The reduced asymptotic front, 
i.e. Afc^{k,Y) x (fcVQsC^^))^" (see Eq. Q), is the straight 
line. 

3. Physical interpretation of universality 

Let us discuss the physical picture underlying the uni- 
versality properties of a general scattering process in high 
energy QCD, beyond the specific example of the BK 
equation. The linear kernels of the saturation equations 
correspond to a developing tree of partons. The tree 
structure is generated by successive branching of partons 
when the rapidity Y increases, which stems from the non- 
abelian character of QCD. For example, the BFKL ker- 
nel is known to correspond to a specific tree of gluons or, 
equivalently at large iVc, of cascading dipoles 19]. 

In the course of the parton branching, the partons get 
more densely packed together, and thus the color field 
gets stronger. The strengh of the field can be qualita- 
tively characterized by the average transverse distance 
d{Y) between neighboring partons. 

The scattering process corresponds to the interaction 
of a source with the partons. As long as d{Y) is larger 
than the resolution 1/Q of the probe, the field seen by the 
probe is weak and thus the interaction merely "counts" 
the partons. This is the dilute regime, in which the evo- 
lution of the observable follows the exponential growth 
of the parton density with Y given by the BFKL equa- 
tion. When d{Y) reaches 1/Q, the color field seen by the 



probe becomes strong, and the probe can interact with 
several partons simultaneously, see Fig. 3. This is the 
saturation transition, and d{Y) plays the role of the sat- 
uration size l/Qs{Y).^ As d{Y) gets even smaller, the 
partons inside the tree may interact and recombine. All 
these finite density or saturation effects generate nonlin- 
ear damping terms in the evolution equations that slow 
down the evolution. Depending on the physical situation, 
the evolution of the interaction amplitude in this regime 
is e.g. described by the BK equation, or by the Balit- 
skii [51 and JIMWLK ^ equations which include more 
saturation effects. 

The structure described here complies with the three 
abovementioned conditions (^jj, (ii), ('ziij for an evolution 
equation to admit traveling wave solutions. The previous 
discussion is general and applies to various modifications 
of the saturation equations. Therefore, we expect the 
universal results presented here to be valid in general, up 
to the replacement of the relevant parameters related to 
the basic tree and which characterize the linear evolution 
kernel. 

4. Conclusion and discussion 

We have pointed out that the structure of the QCD 
parton model is related to general cascading processes 
possessing specific mathematical properties. We have 
obtained all universal terms including new ones in the 
asymptotic expansion of the saturation scale and of the 
unintegrated gluon distribution, see Eqs. H3I4|I . Univer- 
sality means independence of the solution with respect to 
the initial conditions after a long enough rapidity evolu- 
tion and dependence of the solution on a few relevant pa- 
rameters obtained from the kernel of the Zmear evolution: 
7c, xilc), x"(7c) for the saturation scale, plus x'^^Hlc) 
for the front. The results have been obtained from the 
BK equation, but are straightforwardly applicable to a 
whole class of saturation equations. The saturation scale 
and the gluon density above the saturation scale keep 
exactly the same functional form for the universal terms. 

Finally, it is interesting to look for the concrete im- 
prints of saturation on the universal parts of the solution, 
Eqs. (|3I4|1 . They appear in two places: the coefficient of 
the second term of the expansion of the saturation scale 
(| instead of ^ for BFKL, see Eq. Q), and the logarith- 
mic factor in A/'oo in the expression for the asymptotic 
front (Eq. Q), that is responsible for the developing lin- 
ear shape of the reduced front profile in the leftmost part 
of the leading edge, see Fig. 2. This factor appears to be 
a general dynamical property of absorbing walls, near 



This phenomenon is conventionally interpreted as the onset of 
the black disk due to the filling of the phase space by partons of 
size 1/Q 0. 
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FIG. 3: The tree structure. The parton branching is repre- 
sented along the rapidity axis (upper part) and in the trans- 
verse coordinate space at two different rapidities (lower part). 
The probe is represented by a shaded disk of size 1/Q. At ra- 
pidity Yi, the probe counts the partons (linear regime), at Y2, 
the probe counts groups of partons (transition to saturation). 
This effect, and also other effects such as recombination of 
partons inside the tree (r.h.s. of the plot), generate nonlinear 
damping terms in the evolution equations. 



which a linear gradient develops at long times [lH . The 
nonlinearities present in the saturation equations act ef- 
fectively as an absorbing wall in the results of Eqs. H3I4|I . 
This a posteriori justifies the treatment of saturation pro- 
posed in Ref. [ij- 

Numerical solutions of the BK US EH and JIMWLK 
[2H equations are available. The traveling wave behav- 
ior of the solution was seen, and the consistency with 
the first two universal terms in Eq. @ was checked. It 
would be interesting to include the newly found term in 
the comparison between the numerical and analytical so- 
lutions. However, the expansion (|3I4|) is an asymptotic 
series, which is only slowly converging, see Figs. 1,2. 

The question of the phenomenology of saturation is 
still open. For realistic phenomeno logy , one should go 
beyond leading order BK (see e.g. [221) since, as well- 
known, the saturation scale predicted by the leading or- 
der BK equation has a too steep rapidity evolution. 

In a forthcoming publication [23 , we take into account 
some of these subleading effects by extending the dis- 
cussion of universality to the case of the BK equation 



with running coupling. The latter falls in a different 
class of nonlinear equations, but it can be investigated 
through the same general methods, showing the power of 
the mathematical approach. 

We thank Edmond lancu and Dionysis Triantafyl- 
lopoulos for their comments. 
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